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A hybrid version the deformed nuclear potential is suggested, which combines a spherical Woods 
Saxon potential with a deformed Nilsson potential. It removes the problems of the conventional Nils- 
son potential in the mass 130 region. Based on the hybrid potential, tilted axis cranking calculations 
are carried out for the magnetic dipole band in 128 Ba. 



I. INTRODUCTION 

The transitional nuclei in the region A = 130 — 140 show regular AI — 1 regular bands, charac- 
terized by large B(M 1)/B(E2) ratios, the lack of signature splitting and relatively low dynamical 
moments of inertia These bands have an intermediate character. They stand between the collec- 
tive high- if bands of well deformed nuclei, for which collective rotation is the dominant mechanism 
of generating the angular momentum , and the magnetic rotation of near spherical nuclei, for which 
few high-j particles and holes generate most of the angular momentum by means of the shears 
mechanism (see, for example |],^])- The question of how changes magnetic into collective rotation 
has not been studied yet. So far, only the magnetic dipole band in 128 Ba has been investigated 
Q from this point of view. Another intriguing question is the possibility of a chiral character of 
rotation ||. Their softness with respect to triaxial deformations makes the nuclei in the A = 130 
region particularly good candidates for identifying this new symmetry type. 

The Tilted Axis Cranking (TAC) model Q has turned out to be an appropriate theoretical tool 
for the description of the magnetic dipole bands. This model is a natural generalization of the 
cranking model |J for situations where the axis of rotation does not coincide with a principal axis 
of the density distribution of the rotating nucleus, and thus the signature is not a good quantum 
number. Since introduced, TAC has proven to be a reliable approximation for the energies and 
intraband transitions in both normally and weakly deformed nuclei ||. However, in the case of 
the four quasiparticle magnetic dipole band in 128 Ba, the TAC calculations predicted the wrong 
parity and a too early termination of the band and provided only fair description of the electro 
magnetic transition data |p|-p4|. 128 Ba is one of the best studied nuclei in this mass range and the 
above mentioned AI — 1 rotational band is a good test case for the TAC model. 

The purpose of the present work is to try to identify the origin of the discrepancies and remove 
them. The version of TAC used in [Jt| was based on the Nilsson Hamiltonian. The parameters of this 
well-tried potential were carefully adjusted for various mass regions, where it was successfully used in 
standard cranking calculations Jl^] as well as in TAC calculations . However, the A — 130 region 
is known to be problematic for the model. In fact no really satisfactory set of Nilsson parameters 
for this region is available so far ]13| . Thus we attribute the discrepancies of the TAC calculations 
for 128 Ba jjj with the later measurements [p[-p4| to the general problems of the Nilsson potential in 
this region. 

On the other hand, the Woods-Saxon potential works very well around A — 130 . Encouraged 
by this, we adapt the Nilsson potential as close as possible to the Woods-Saxon one. We call this the 
hybrid potential, which is the basis for new TAC calculations for 128 Ba. Instead of the parameterizing 
the single particle levels of the spherical modified oscillator in the standard way by means of an Is- 
and an Z 2 -term, the hybrid model directly takes the energies of the spherical Woods-Saxon potential. 
The deformed part of the hybrid potential is an anisotropic harmonic oscillator. This compromise 
keeps the simplicity of the Nilsson potential, because coupling between the oscillator shell can be 
approximately taken into account by means of stretched coordinates Jl^ , and it amounts to a minor 
modification of the existing TAC code. On the other hand, it has turned out to be a quite good 
approximation of the realistic flat bottom potential as long as the deformation is moderate. The 
hybrid was used to calculate the triaxial shapes of liquid sodium clusters jlf|. The results agree 
very well with later calculations using the correct radial profile of the deformed part of the potential 
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II. TAC IMPLEMENTATION 



The TAC model is discussed in more detail in [|2|-|5|. The brief presentation in this paper focuses 
at the suggested improvement of this approach. The starting point is the mean field Routhian[] 

h' = h Bph + V def (e 2 ,i, £i) " A (P + P + ) - XN - lu (sintfji + costf j 3 ) (1) 

where h sp u denotes the spherical part including the spin-orbit term and Vdef{£2, 7, £4) the deformed 
part of the Nilsson single-particle Hamiltonian. (see, e.g. (l^]). The pairing field in eq.(l) is deter- 
mined by the gap parameter A and the monopole pairing operator P = CkCj: while the chemical 
potential XN is needed to satisfy on average the particle number conservation. The two-dimensional 
cranking term u> (sin ftji + cos^j'3) is the new element of the TAC, as compared to the standard 
cranking model, which is recovered for ft = or 90° . The angle ft fixes the tilt of the cranking 
axis with respect to the intrinsic 3-axis within the principal (1-3) plane of the deformed potential. 
By diagonalization in stretched coordinates, neglecting AJV = 2 shell mixing, this TAC Routhian 
yields quasi-particle energies and quasi-particle states, from which the many-body configuration 
\uj, £2, 7, £4, ft > of interest is constructed. 

The mean field is found for a given frequency lu and fixed configuration by minimizing the total 
Routhian 

E'(u, £2, 7, £4, ft) =< w, £2, 7: £4, ft I h' lu, £ 2 , 7, £4, ft > (2) 

with respect to the deformation parameters (£2, 7, £4) and the tilt angle ft. The value of A is kept 
fixed at two values: 80% of the experimental odd-even mass difference and zero. At the equilibrium 
angle ft = ft (minimum) the cranking axis is parallel to the direction of the angular momentum 
vector J = (< ji >, < j3 >). After the minimum is found the various electro magnetic observables 
of interest are obtained by means of the following semiclassical expressions 

< I — 21 — 2\M-2(E2)\II >=< -M_2(£2) >= (3) 
= Vsf(!t) [v 7 ! <Q'o > (sin ft) 2 + j <Q' 2 + Q'-2 > (l + (costf) 2 )] , (4) 

< / - II - l\M-i(E2)\II >=< M-i{E2) >= (5) 

= V / 3?(f) [ Sil1 ftcOsft(y/¥ 2 <Q' Q > -\<Q'2 + Q'-2>)\ ■ (6) 

<I -II- l\M-i(Ml)\II >=< M-i(Ml) >= (7) 

= ViF I^ 3 sin ^ ~ a* 1 cos ^] ' ( 8 ) 

where < M v > is the expectation value of the transition operator with the TAC configuration 
I >, the components of which refer to the lab system. The intrinsic quadrupole moments Q'^ are 
calculated with respect to the principle axes (1, 2, 3). The same holds for the magnetic moments 

Mi = /J-n(Ji, p + (775.58 - 1) Si, p - ?73.82 Si,„), fj,3 = A*iv(J3,p + (775.58 — 1) Ss, v — ??3.82 S^n), (9) 

where the free nucleonic magnetic moments are attenuated by a factor of 77 = 0.7. The reduced 
transition probabilities are 

B(M1, AI = 1) =< M-x{Ml) >' 2 
B(E2,AI = 2) =< M-2(E2) > 2 

The mixing ratio is 

<M- l (E2)> 

<M-x{Ml)> K ' 

We apply the Strutinsky renormalization procedure to calculate the total Routhian 

E'(w) = E ld (lu = 0) - E smooth + <(j\h'\w>, (11) 



x For simplicity, only one type of particles is spelled out. The extension to both types is obvious. 
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where Eld = Eld{£2,"/, £4, £4) means the liquid drop energy and E amoot h is the the smooth part of 
the mean-field energy calculated from the single-particle energies at ui = 0. This version of the TAC 
has turned out to be quite successful for well deformed nuclei (see j^-^] and references therein). 

The new element of the present paper consists in the hybrid potential, which approximates the 
well-established deformed Woods-Saxon potential, yet preserving the existing convenient TAC envi- 
ronment. For this purpose the spherical part h sp h in Eq.(l) is replaced the spherical Woods-Saxon 
Hamiltonian for the nucleus of interest. In the present work the universal Woods-Saxon parameters 
are used (see e.g. p7|). 

Technically, the replacement is rather simple, because the existing TAC code uses states of good 
l,j,m as a basis. The spherical Nilsson energies ■ axe replaced by the spherical Woods - Saxon 

energies e^f' '•. It turns out to be unproblematic to associate the quantum numbers of the two differ- 
ent potentials. For a given combination l,j the third quantum number N is found by counting from 
the state with the lowest energy. The fact that the spherical Woods-Saxon code uses a harmonic 
oscillator basis permitted a check of the algorithm. The major component of the Woods-Saxon 
wavefunction agrees with the state found by our counting algorithm. In the high-lying part of the 
single-particle spectrum (three shells above the valence shell or higher) there are occasional ambi- 
guities in assigning the states. Small errors of this kind are not expected to have any consequences 
at moderate or small deformation. The states do not couple strongly to the states near the Fermi 
surface. They contribute only to the smooth level density used in the Strutinsky renormalization, 
which will not be affected by small shifts of the levels. 

Such a replacement of the spherical single-particle energies is a common practice in large scale 
shell-model configuration mixing and similar calculations, where they are often used as adjustable 
parameters The effect of the replacement is illustrated on Fig.l, which shows the deformation 
dependence of the proton single-particle levels of 128 Ba for the Nilsson, the Woods-Saxon and the 
hybrid Hamiltonians. The close similarity of the levels of the Woods-Saxon and the hybrid models 
is obvious. The hybrid has a somewhat later and sharper crossing between the positive parity levels, 
which also show stronger tendency to arrange into pairs of pseudo spin doublets. These treats are 
inherited from the Nilsson Hamiltonian, which controls the change with deformation. The main 
difference between the Nilsson Hamiltonian and the other two is the lower energy of the negative 
parity levels originating from /in/2- It seems to be the reason for the discrepancies between the 
previous TAC predictions and the experiment, as will be demonstrated in the next section. 



III. THE Ml BAND IN 12S BA 

In the previous TAC calculations for the A7 = 1 band in 128 Ba Q the four quasi-particle con- 
figurations [n(h 11 / 2 ) 2 i / {hii/2(d 5 /2g7/2))] f^] for the negative parity and [n(h 11 /2) 2 v(h 11 / 2 ) 2 ] for the 
positive parity were found to be the lowest ones in energy at the oblate deformation of e = 0.26, 
7 = 60°. This deformation was determined at u) = 0.2 MeV by minimizing the total Routhian 
calculated from the Quadrupole-Quadrupole interaction. With the present version of TAC we find 
a significantly smaller prolate equilibrium deformation of e = 0.205 and 7 = 0°. The lowest four 
quasi-particle configuration now turns out to be [7r(/in/2(d5/237/2))i / (/iii/2(rf5/257/2))]- Fig. 2 shows 
the quasi proton and quasi neutron levels for various cranking frequencies at tilt angle ■& — 90° and 
for various angles at u) — 0.2 MeV. The equilibrium value of #, which minimizes E' , is found to be 
1? = 52.5° for u = 0.2 MeV. 

As seen, the different position of the orbitals in the hybrid TAC has drastic consequences. 

The deformation changes from oblate to prolate, resulting in a different configuration of the Ml 
band. This is not surprising in a region where the energy difference between oblate and prolate 
shape is small. 

The calculations of the AI = 1 band in the present work are built on this new configuration 
[7r(/in/2(rf5/2<77/2))(j / 'iii/2(d5/2<77/2))], which is the lowest four-quasi particle TAC solution. In agree- 
ment with the experiment, it has positive parity. The excitation of four quasi-particles significantly 
reduces the pairing gaps. In order to better grasp the influence of this blocking on characteristics 
of the band, we did two calculations: one with A„ = 0.88 MeV and A,r = 1.04 MeV corresponding 
to 80% of the experimental even odd mass difference and one with zero pairing. 



indicating the major components, we denote the mixed Nilsson state by ((^5/237/2)- 
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Electro magnetic transition properties present a stringent test of the nuclear models. The 
experimental information about lifetimes and mixing ratios of the magnetic dipole band was 
accumulated in several experiments |fl-alj|. Fig. 3 compares the experimental B(M1) and 
B(E2) values and mixing ratios with the results of TAC model obtained for the configuration 
[7r(/i 11 / 2 (rf5/2<?7/2)) I/ (' 1 ii/2(d5/2<?7/2))]- All electro-magnetic characteristics of the band are well re- 
produced by both the paired and unpaired calculations. Curiously, the quenching of pairing influ- 
ences to some extent the B(M1) values but leaves almost unchanged the B(E2) values up to spin 
18h. The TAC calculation seems to slightly overestimate the deformation. 

In contrast with the previous TAC calculation, it is possible to follow the band all the way up 
to spin 26h. Fig.4 shows the measured and the calculated function J(lj) of the spin on the angular 
frequency. It is more sensitive to the changes of the pair correlations. While the unpaired calculation 
gives a nearly linear function with the moment of inertia = dj/duj close to the measured one, 
the paired calculation exhibits a substantially lower moment of inertia for low rotational frequencies 
and an upbend at higher ones. The experiment is in between. This seems to indicate that the pair 
field is weak in this nucleus and a more refined treatment of pairing is needed. It is noted that in 
the calculations the band extends down to u> — 0.1 MeV and J = llh. In experiment there is an 
irregularity around I = 12, 13h. It may be caused by mixing of the 13 + state with the another 13 + 
state, which lies nearby in energy (AE = 0.045 MeV) and into which the 14 + also decays jll| . 

In the A J = 1 band was analyzed in terms of a pure high-K band using the familiar 

expressions for the B(Ml) and B(E2) values for the axial symmetric rotor [Q. Adjusting three 
free parameters, the K- value, the intrinsic quadrupole moment Q , and the gyromagnetic factor 
Iff A" — gii\, a good fit of the electro magnetic decay data was obtained. The quality is practically 
the same as in our calculation without parameters in Fig.4. The TAC calculation contains much 
more physical information as e.g. the specific configuration on which the band is built and the band 
energies. The knowledge of the intrinsic state can be used to derive further structure information, 
e. g. the geometrical coupling scheme shown in Fig. 5, which enables one to see how the total 
spin is formed from the quasiparticle orbitals and how it changes with the rotational frequency. 
Apparently, most of the angular momentum gain along the band is of collective nature, however 
the high-j quasiparticles from proton and neutron (££5/257/2) an d /111/2 orbitals do also substantially 
contribute by means of the shears mechanism. With increasing frequency, the 3-component of the 
spin vector J stays practically at < J3 >w 9h. However, this does not mean that the corresponding 
TAC configuration behaves like a structureless the high-K rotor. Fig. 6 shows that the calculated 
dependence of the tilt angle on the rotational frequency by no means follows curve expected 

for the strong coupling limit. Thus, the present case lies in between a good shears band and a good 
high-K band. The A7 = 1 band in 128 Ba is an example for a rotational band of intermediate nature. 



IV. SUMMARY 



The strength of the TAC model is that it can predict the appearance of A7 = 1 rotational bands 
and is able to describe microscopically their electro magnetic decay properties. This is achieved 
by taking into consideration the orientation of the rotational axis with respect to the deformed 
potential, which is fixed along a principal axis in the conventional cranking model. The intrinsic 
TAC configuration of a rotational band is found by searching for a local minimum on the multi- 
parameter surface of the total Routhian. Therefore, it is crucial to calculate these surfaces as 
reliable as possible. In the present work the Strutinsky renormalization and a hybrid single-particle 
potential were implemented in order to improve the total Routhian. The hybrid potential combines 
the spherical Woods-Saxon single-particle energies with the deformed part of the Nilsson potential. 
For moderately deformed nuclei it gives deformed single-particle levels that are quite close to the 
deformed Woods-Saxon levels. 

We applied the hybrid TAC to the previously investigated four quasiparticle magnetic dipole band 
in 128 Ba. The lowest equilibrium configuration is found to be [7r(/iii/ 2 (d5/257/2))j / (^ii/2(rf5/2fl7/2))]- 
It has positive parity and a prolate axial deformation of £2 = 0.20. The microscopic TAC calculations 
describes rather well the experimental energies, B(M1) and B(E2) values, as well as the branching 
and mixing ratios. The dipole band in 128 Ba has an intermediate structure. A comparable amount 
of angular momentum is generated by the shears mechanism, active for the /in/2 and (1^5/2(77/2) 
quasiparticles, and by collective rotation. It is an example for the transition from collective to 
magnetic rotation. 

The hybrid potential turned out to be crucial for the good agreement between the calculation and 
the data. It substantially improves the Nilsson potential, which has problems in the region around 
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mass 130. It seems to be a promising starting point for studying the intriguing interplay between 
triaxial deformation and the orientation of the rotational axis. 
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FIG. 1. Proton single-particle levels as function of deformation for the Woods-Saxon (bottom), hybrid (middle) and Nilsson 
(top) model, respectively. Full lines present positive parity states, broken lines - negative parity states. The Fermi energy lies 
at about 51.3 MeV for the Woods-Saxon and the hybrid models and at 45.6 MeV for the Nilsson model. 
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FIG. 2. Proton (a) and neutron (b) quasi-particle levels as function of rotational frequency Too and tilt angle Full 
lines correspond to positive parity, broken lines - to negative parity states. In both cases the lowest two quasi-particle 
configurations involve one predominantly (^5/2^7/2) an d one predominantly /111/2 state. The deformation parameters are 
£2 = 0.205, £4 — —0.01, 7 = 0°, which are the equilibrium values at ui = 0.2 MeV. 
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FIG. 3. Electromagnetic observables as obtained by paired calculation (full lines), unpaired calculation (dotted lines) and 
experiment. 
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FIG. 4. Dependence of the angular momentum J = / on the rotational frequency hco as obtained by paired calculation (full 
line) , unpaired calculation (dotted line) and experiment (crosses) . The experimental frequency is extracted from the measured 
7— energy data using the relations hu(I) — — E(I) — E(I — 1). 
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FIG. 5. Spin decomposition in terms of individual quasi-particle contributions for two different rotational frequencies. 
Quasi-particle configurations are denoted by their predominant components. 
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FIG. 6. Equilibrium tilt angle ■& as a function of the rotational frequency Tiuj. Full line - paired calculation; dotted line - 
unpaired calculation; broken-dotted line - strong coupling limit result. 
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